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We discuss the spin-wave interaction in two-dimensional (2D) Heisenberg ferromagnet (FM) with 
dipolar forces at Tc 2> T > using 1/S expansion. A comprehensive analysis is carried out of 
the first 1/S corrections to the spin- wave spectrum. In particular, similar to 3D FM discussed in 
our previous paper A.V. Syromyatnikov, PRB 74, 014435 (2006), we obtain that the spin- wave 
interaction leads to the gap in the spectrum tk renormalizing greatly the bare gapless spectrum at 
small momenta k. Expressions for the spin-wave damping Fk are derived self-consistently and it is 
concluded that magnons are well-defined quasi-particles in both quantum and classical 2D FMs at 
small T. We observe thermal enhancement of both Fk and Fk/tk at small momenta. In particular, 
a peak appears in Fk and Fk/tk at small k and at any given direction of k. If S' ~ 1 the height of 
the peak in Fk/ck is not larger than a value proportional to T/D <C 1, where D is the spin-wave 
stiffness. In the case of large spins S ^ 1 the peak in Fk/tk cannot be greater than that of the 
classical 2D FM found at = which height is small only numerically: Fo/eo ~ 0.16 for the simple 
square lattice. Frustrating next-nearest- neighbor exchange coupling increases Fo/eo in classical 2D 
FM only slightly. We find expressions for spin Green's functions and the magnetization. The latter 
differs from the well-known resuh by S.V. Maleev, Sov. Phys. JETP 43, 1240 (1976). The effect 
of the exchange anisotropy is also discussed briefiy. Higher order corrections to the spectrum are 
considered and it is concluded that they are small compared to the first corrections obtained. Our 
results contradict to findings of the previous works (Ar. Abanov, A. Kashuba, V.L. Pokrovsky, PRL 
77, 2554 (1996); PRB 56, 3181 (1997)) in which a diffusion spin-wave mode was obtained at small 
momenta. It is shown that the origin of this discrepancy is that the spin-wave gap was ignored in 
the previous studies. 

PACS numbers: 75.70.Ak, 75.30.Ds, 75.10.Jm, 75.10.Dg 

I. INTRODUCTION 

Magnetic properties of thin films and (quasi-) two-dimensional magnetic materials are of great interest nowii'^ This 
interest is stimulated by recent advances in film growth techniques and numerous technological applications of magnetic 
films, including uses in electronics and data storage. A realistic theoretical model of low-dimensional magnetic systems 
must include the exchange interaction, the dipolar interaction, and the magnetocrystalline anisotropy.A Despite its 
smallness the long-range dipolar interaction plays essential role in 2D magnets. In particular, it violates the Mermin- 
Wagner theorem'^ and leads to stabilization of the long-range magnetic order at finite temperature in 2D magnets^^J^i^ 
Some peculiar features have been observed recently at T ^ related with the dipolar interaction both in 2D and 
3D Heisenberg ferromagnets (FMs). 

In 3D FM the problem of infrared singularities arose. It was obtained in Ref.^ that dipolar forces lead to strong 
long-wavelength fluctuations manifesting themselves in infrared divergence of the first perturbation corrections to the 
uniform longitudinal spin susceptibility: x\\i^ ^ 0) ~ iT/uj. An infrared divergent contribution to the spin-wave 
stiffness was obtained in Ref.^ as a result of analysis of the first 1/5 corrections to the spin- wave spectrum. Thus, 
the problem arose of analysis of the whole perturbation series in order to find spin-wave spectrum and longitudinal 
spin susceptibility at small momenta. Appearance of the infrared singularities in these papers is related to the fact 
that the spectrum is gapless in 3D FM in the spin-wave approximation.® 

First perturbation corrections to the spin-wave spectrum were analyzed in classical 2D FM with dipolar forces in 
Refs. It was found that the imaginary part of these corrections exceeds the bare gapless spectrum at small enough 
momenta k. Then as a result of self-consistent calculations a diffusion spin-wave mode was obtained at very small k. 
It was argued in Refs.^'^° that despite the analysis is carried out for classical 2D FM the diffusion mode should be 
observed also in quantum 2D FM. 

Meantime we find that quite an unusual property of ferromagnets with dipolar forces is ignored in the previous 
studies that is crucial for the reported peculiarities — interaction between magnons leads to appearance of a gap in 
the spin-wave spectrum that renormalizes the bare gapless spectrum greatly at small momenta. Thus, we obtain such 
gap in 3D FM in the first order of 1/S in our recent paper'^l . It was found to be proportional to wq^/ Suoq/ J sin 0k, 
where 



(1) 
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is the characteristic dipolar energy, vq is the unit cell volume, is the angle between momentum k and magnetization 
and J is the exchange value. We show that this gap screens the infrared singularities obtained in RefSiSiZ and first 
perturbation corrections to the observables found self-consistently are small. Naturally, one can expect existence of 
such gap in 2D FM too. 

It should be noted that appearance of the gap in 2D and 3D FMs is quite expected because dipolar interaction due 
to its symmetry and long-range nature violates the Goldstone theorem. Besides, the spin-wave gap was observed in 2D 
antiferromagnet with dipolar interaction in the zeroth order of 1/5* (i.e., in the spin-wave approximation).^ Existence 
of the gap in the spectrum of 3D FM was anticipated long times ago. It is well known that within the first order of 
1/S dipolar and pseudodipolar forces lead to anisotropic corrections to the total energy of 3D FMji^ii^iii As a result, 
directions along edges of the cube are energetically favorable in a simple cubic lattice whereas the magnetization 
should be parallel to a body diagonal of the cube in the face-centered cubic lattice and the body-centered cubic 
latticeJ^ii^'i^ii^ It has been pointed out by Keffer^'''^^ that the anisotropic terms in the total energy of a ferromagnet 
should be accompanied with an "energy shift" in the spin- wave spectrum. We confirmed this long-standing statement 
for 3D FM in our previous paper— and demonstrated the relation between the gap and the anisotropy at T = 0. A 
four-fold in-plane anisotropy in square 2D FM caused by quantum fluctuations was also observed before and the value 
of this anisotropy was investigated numerically in Ref.^^ for T ^ SJ. 

The situation is slightly different in classical FMs with dipolar interaction. In particular, the ground state remains 
infinitely degenerate according to the rotation about the axis perpendicular to the plane in the classical 2D FM. 
But it is well known that in classical systems with degenerate ground states some of these states can be selected via 
order-by-disorder mechanism at finite temperature fiii^iii^ Order-by-disorder effect was demonstrated in 2D systems 
of classical spins confined to lie within the plane and coupled via (i) short-range dipolar-like interaction (Ref?i2) and 
(ii) long-range dipolar interaction (Refsi^i^i). Then, an in-plane anisotropy arises at finite temperature. Notably, 
it was found in Ref.^^ that the thermal selection of the ground state is accompanied by the appearance of the gap 
proportional to T in the spin- wave spectrum originating from the spin-wave interaction and leading to the finite value 
of the order parameter. Then, one expects also appearance of the spin-wave gap in classical 2D FM with dipolar 
forces. 

In the present paper we carry out a comprehensive analysis of the first 1/S corrections to the spin- wave spectrum 
in 2D FM with dipolar interactio n. Simi lar to 3D FM we ob tain the spin-wa ve gap A in the spectrum Ck which 
appears to be proportional to uoq^Sujq/J for S" ~ 1 and SujQ^/{uJn/ J) (T/S'^J) for S" > 1 and Tc » T > SJ. This 
gap renormalizes greatly the bare gapless spectrum at momenta k < {A/ (Sujq))'^ . The limiting case of classical spins 
is also discussed at T ^ j and the gap was found to be proportional to w^Tw j p ^ where j and w are values of the 
exchange and the characteristic dipolar energy in the classical model, respectively. It is demonstrated below that the 
spin-wave gap that was ignored in Refs.^'^" are much larger than the energy of the diffusion mode obtained in those 
papers. We derive the spin- wave damping Fk self-consistently and find that spin waves are well-defined quasi-particles 
in both quantum and classical 2D FMs at small T. 

Interestingly, we observe thermal enhancement of the damping at small momenta. In particular, in quantum 2D 
FM we obtain a peak in both Fk and Fk/ck at fc <C wq/ J, T ^ Sloq and at any given k direction. If S" 1 the height 
of the peak in Fk/ck cannot axceed a value proportional to T/{SJ) <C 1. In the case of large spins 5* 3> 1 the peak 
in Fk/ek cannot be greater than that of the classical 2D FM found at /c = which height is small only numerically: 
To/fo ~ 0.16 for the simple square lattice. Frustrating next-nearest-neighbor exchange coupling increases Fo/eo fi^ 
classical 2D FM only slightly. 

We derive expressions for spin Green's functions and the magnetization. The latter differs from the well-known 
result of Ref4. The effect of the exchange anisotropy is also discussed briefly. Higher order corrections to the spectrum 
are considered and it is concluded that they are small compared to the first corrections obtained. We derive the four- 
fold in-plain anisotropy in the total energy of quantum 2D FM that makes directions of the magnetization along edges 
of the square to be energetically favorable in the simple square lattice. We also demonstrate the relation between this 
anisotropy and the spin-wave gap at T = as it was done in our previous paper^^ for 3D FM. 

The rest of the present paper is organized as follows. The Hamiltonian transformation and the technique are 
discussed in Sec. |lll First 1/ S corrections to the real and imaginary parts of the spin- wave spectrum are considered 
in Sees. [IlT] and |TVl respectively. In Sec. |V] we study the case of large spin values and consider the limit of classical 
spins. In Sec. IVII we i) demonstrate the relation between the anisotropic term in the total energy appearing due to 
the dipolar interaction and the spin-wave gap at T = 0, (ii) discuss the further order \/ S corrections to the spin- wave 
spectrum, (iii) calculate the magnetization taking into account the spin-wave spectrum renormalization, (iv) derive 
the spin Green's functions in the first order of 1/S, and (v) discuss briefly the effect of the easy-plane anisotropy and 
consider the spectrum renormalization obtained in Refs.^"-^". Sec. IVIll contains a detailed summary and our conclusion. 
Three appendixes are included with some details of calculations. 
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II. HAMILTONIAN TRANSFORMATION AND TECHNIQUE 



The Haniiltonian of a ferromagnet with dipolar interaction has the form 



n 



yl3 



2 ^ 
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Taking the Fourier transformation we have from Eq. 

1 ^ 
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where Jk = J2i Jim exp(ikR/„i) and = J2i Qim exp(?kR;m)- We direct y-axis perpendicular to the lattice as is 
shown in Fig. [T] Dipolar tensor possesses the well-known propertiesi'"* at fc ^ 1 , which are independent of the 
lattice type and the orientation of x and z axes relative to the lattice, 



Ql^ = Wo ( -^a + '- 1 Syi3, where (3 = x,y, z, 



2 fc2 



, where p,(3 = x, z, 



a = 



1 



(5) 
(6) 
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where we set the lattice spacing to be equal to unity, ujq is the characteristic dipolar energy given by Eq. ([T]) and a 
is a constant that is equal approximately to 1.078 for the simple square lattice. It is seen from Eqs. ^ and ^ 
that dipolar forces lead to easy-plane anisotropy in the energy of the classical 2D FM with y to be a hard axisji We 
show in the next section that quantum and thermal fluctuations lead also to an in-plain anisotropy. We direct z-axis 
along the uniform magnetization as is shown in Fig. [T] 
After Dyson-Maleev transformation 
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Hamiltonian (j4|) has the form Ti = Eo + Yl^i=i '^i^ where Eq is the ground state energy and Hi denote terms containing 
products of i operators a and . One should take into account terms up to 7^4 to calculate corrections of the first 
order in 1/5. Tli — because it contains only Qq^ with p ^ j3. For the remaining necessary terms one has 
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where we drop index k in Eqs. (fTO]) and (|lip . is the number of spins in the lattice and 
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where D is the spin-wave stiffness, (/)k is the angle between k and the magnetization and the expressions after are 
approximate values of the corresponding quantities at A: <C 1. One has for the coupling between only nearest neighbor 
spins on the simple square lattice D ~ SJ. In the spin-wave approximation wc find for the magnon spectrum 



= /b^^ ^(Dk^ + ^^o«) (^Dk^ + sin^ 0k ) (14) 

in accordance with the well-known result.^ 

To perform the calculations it is convenient to introduce the following retarded Green's functions: G(ciJ,k) — 
{ak,al}uj, F{uj,'k) = (ak,a-k)(^, G(cj,k) = (ali^,a_k)w = G*(-a;,-k) and F\uj,k) = (al,^., aj^)i^ = F*{-uj,~'k). 
We have two sets of Dyson equations for them. One of these sets has the form 

G{Lu,k) = G(")(tJ,k) -|-GW(tJ,k)E(w,k)G(a;,k) +G(°)(^,k)[Bk + n(cj,k)]Ft(w,k), 

F^iuj, k) = G*°^ {uj, k)S(c^, k)^^t (^^ k) + g'"' {lj, k) [Bk + nt (c^, k)]G{uj, k), ^^^^ 

where G'"^ (w, k) = (tj — i?k + j(5)^^ is the bare Green's function and E, E, 11 and H^' are the self-energy parts. Solving 
Eqs. ([15]) one obtains 



G{LU,k) 

F(u;,k) 



! + E^ + E{Lu,k) 

Vicj,k) 
gk + n(co,k) 
2?(^,k) ' 



G(w,k = — — — , 16 

V{uj,k) 



where 



P(w,k) = (w-f i(5)^-e^-0(w,k), (17) 
n{uj,k) = i;k(i; + E)-Bk(n-hrf)-(u;-|-i(5)(E-E)-ntf+EE, (18) 

and ek is given by Eq. ((M]). Quantity il.{u!,k) given by Eq. (fTS]) describes renormalization of the spin- wave spectrum 
square. We calculate the real part of 17 (w, k) in the next section and analyze its imaginary part in Sec. IIVI The last 
two terms in Eq. p8|) give corrections of at least second order in 1/5* and arc not considered in Sees. IIIIfTvl We imply 
that S' ~ 1 in the next two sections and discuss large S in Sec. |Vl 



III. RENORMALIZATION OF THE REAL PART OF THE SPIN- WAVE SPECTRUM 

Corrections to the spin-wave spectrum to be obtained are proportional to sums over momenta in which summands 
depend on the components of the dipolar tensor ■ In some of these sums summation over small momenta is 
important and one can use expressions © and ^ for Q^^ ■ Meantime there are sums in which summation over large 
momenta is essential and which, consequently, depend on the direction of the quantized axis within the plane and the 
lattice type. Thus, one should bear in mind what is the direction of magnetization in the ground state. 

Therefore the well-known fact should be taken into account that dipolar and pseudodipolar interactions lead to the 
dependence of the energy of a ferromagnet on the direction of quantized axis,^-- We show now that, being established 
first for 3D FM, this finding remains also valid for 2D FM: apart from the easy-plain anisotropy found in Ref.* and 
discussed above there is also an in-plain anisotropy. As in 3D FM, the first l/S'-correction to the classical energy Eq 
having the form 
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gives rise to such in-plain anisotropy, where the isotropic term is omitted in the right part of Eq. (jl9p . After direct 
calculations one obtains for a square lattice 
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where 7^ are direction cosines of the magnetization and components of the dipolar tensor in Eq. (j2ip are taken relative 
to square axes. The constant C should be calculated numerically because summation over large momenta is important 
in Eq. (pij) and one cannot use Eqs. ([5]) and (O for the dipolar tensor components. This calculation can be carried 
out using the dipolar sums computation technique (see, e.g., Refi^ and references therein) with the result C ~ 0.0082 
for exchange coupling between only nearest neighbor spins on the simple square lattice. Because C > an edge of 
the square is the easy direction. Notice also that the value of the dipolar anisotropy in a cubic 3D FM is proportional 
to a constant that also has the form (|2ip . where, naturally, summation is taken over a 3D lattice. The in-plane 
anisotropy caused by dipolar interaction was investigated numerically aX T ^ D in Ref.^^. 

We study now separately bubble diagrams shown in Fig.[2l^a) and (b), and the loop diagram presented in Fig.^c). 



A. Bubble diagrams 



Let us start with the diagram shown in Fig. [2Ja). It appears from three-magnons terms (|10p and gives zero. To 
demonstrate this we make all possible couplings of two operators a and a) in Eq. (jlOp : 



a\ + flo 
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2^ ^ 2e„ ^1 
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where N^^ — {e'^'^/^ — 1) ^ is Plank's function. Expression (|22p is equal to zero because — — Qqf , -Eq = -Eq' and 
■Bq = -Bq', where q = [q^, qz) and q' = {-qx, qz)- 

The Hartree-Fock diagram presented in Fig. [Hb) comes from 7i4-terms given by Eq. pT|) . After simple calculations 
we obtain for the contribution to Q,{lo, k) from this diagram 
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At zero temperature TVq = in Eq. ()23p . In this case the spectrum is renormalized by quantum fluctuations only. As 
the temperature increases, corrections from terms in Eq. (j23p containing TVq become larger. They exceed terms in 
Eq. ((23| not containing A^q above a certain temperature. We find below that this temperature is of the order of Sloq. 
Then, it is convenient to discuss separately regimes T <ti Sloq and T 3> 5^0- 



1. T^Sujo 

Taking A^'q ~ in Eq. ([23|) . we obtain in the leading order of 



XX 

q 



201 



(24) 



Here the first term is of the order of fc^Wp ln(Z?/ [Sujq)). Then it gives positive negligibly small correction to the bare 
spectrum (jl4p . In contrast, the second term in Eq. (|24p . being independent of k, contributes to the spin-wave gap. It 
is much greater than the bare spectrum at /c <C Sloq/D. 



6 



2. T > St^o 



Terms in Eq. (|23p containing A^q come into play at such T and we have in the leading order of cjq 



n('')(cj,k) 



{Dk^f '^WiT) - Dk'^uj(,a[W{T) + V{T)] 



^0 — t/a 



Jo 
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where w — (IGtt) ^ dkk/{e'' — 1) = 7r/96 and w = l/(27r). The first two terms in Eq. (|25p do not change the 
structure of the bare spectrum resulting in the following renormalization of constants D and ujq in Eq. (|14p : 



D 



1 - -WiT) 



U!o ^ LUq 



1 - gV{T) 



(28) 



In contrast, the last term in Eq. ([25)1 changes the form of the spectrum and contributes to the gap. Notice that the 
gap in Eq. ((25| has the same form as in Eq. (|24p . Thermal corrections to the gap are small being of the order of 
u;lT\n{T/{ScJo))/D\^» 

Comparing first two terms in Eq. ([25]) with the first term in Eq. (|24p one infers that thermal fc-dependent corrections 
become much larger than quantum ones at T 3> Stoo- 



B. Loop diagram 

We turn now to the loop diagram shown in Fig. [DJc). It originates from Tia-terms ([TU]) in the Hamiltonian. As a 
result of simple but tedious calculations some details of which are presented in Appendix[X]we have for the contribution 
to the real part of n{uj, k) from this diagram at /c <C 1 

Rel7(3)(^,k) . -^5:1^, (29) 

where we set k = under the sum because summation over large q is essential. Notice that we discard in Eq. (|29p all 
the terms that are much smaller than O^'*^ (uj, k) given by Eqs. (p4|l and ([25]) . Temperature corrections to Rell'^-' (w, k) 
are negligible. 



C. Resulting expressions 

One can derive now the resulting expression for Reri(w,k) using Eqs. ([24| . ([25| and ([29]) . 



1. T<Sajo 



We obtain from Eqs. ^ and ([29l) 



Rel](c.,k) ^ Dk^^Y. ^-\ ^^% A^ (30) 

q ^ 



aCS^"^, (31) 



where C is given by Eq. ((2T|) . As it is mentioned above, C « 0.0082 for the simple square lattice with the exchange 
coupling between nearest neighbor spins only. The first term in Eq. ([30]) originates from ri('')(cj, k) and to the second 

term contribute both n^^^iuj,^.) and f^('')(w,k). The spin- wave gap A given by Eq. (I3ip is proportional to cJq^^. To 
illustrate this result we plot in Fig. [3] renormalized and the bare spin- wave spectrum for = 0, (/)k = 7r/4 and 
'/'k = 7r/2 assuming that there is exchange coupling between nearest neighbor spins only, ujq = 0.05 J and S = 1/2. 
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2. T > Suj^ 

One has from Eqs. ^ and ^ 

Ref](w,k) = -{Dk^)^\w{T)- Dk^ujQO^(T)^V{T)\^- S^, (32) 

where PF(r), F(T) and A are given by Eqs. (|26p . (jST]) and (|3ip . respectively. Notice that the spin-wave gap A has 
the same form as at T ^ StJo, thermal corrections to its square are negligible being of the order of ujqT^/'^ / D^^"^ . 

We infer comparing Eq. with Eqs. ([50)1 and that renormalization of the bare spectrum is small at T ^ £> 
and k 3> Suq/D. In contrast, the spectrum renormalization is significant at smaller k due to the gap that is much 
larger than all fc-dependent terms at fc < {A/{SllJo))'^ . We use renormalized spectrum below for self-consistent 
calculation of the spin- wave damping and estimation of higher order 1/5* corrections. 



IV. SPIN- WAVE DAMPING 



We discuss in this section the imaginary part of r2(a;,k) to which only the loop diagram contributes shown in 
Fig. mjc). Corresponding calculations are rather cumbersome and we discuss only results here. One refers to Ap- 
pendix [X] for some details of the calc ulations . We discuss 2D FM on t he simple square lattice in this section. Three 
regions should be considered: k ^ Sujq/D, Sloq/D <C fc <C ^JSujq/D and k < (A/(S'aJo))^ at which the real part of 
the spectrum has the from ek ~ Dk^ , ~ k^/ceSDuo and ek ~ A, respectively. Im51(ti;,k) is an odd function of uj 
and we calculate it for uj = only. The spin-wave damping Fk at momentum k is found below using the relation 

Imfljuj ek, k) 

rk = . (33) 



fc> ^/Sluq/D 



One obtains for k 3> \J Sujq / D 

Tk = Ek 

where 
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(34) 



(35) 



and A is given by Eq. (j3ip . It is seen that the damping is anisotropic in this regime: it is smaller along edges of the 
square and it reaches maxima along diagonals of the square. 



B. Sluo/D < fc < y^Scoo/D 



One has for Sujq/D < A: < ^JS^o/D 
Tk 
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where 



A k Z"^^"^' dq ^^^^ 

M^) - 4 / " (38) 

J-^T^ (1-9^) V'C((7,k) 

/■"^ a + 1 1 

^3(k) = 2 dq l^ (39) 



15 /-v"!^ g2 



e,(k) = - / dq— -^-2^/:(g,k)(^^-^^/^ z = 1,2,3, (40) 
4 J-v^T^ (l-f? ) 



8. = 4^3/ .,<i±i2!{;#i?.5,31, (41, 



(l + 2g2)4 

£(g,k) - (l-g2)2_^|i_^2|_^(^2^3)^ (42) 



2 



£1 
£2 



- 0' + M) , (45) 
^ (/9 + e + + 0' + 16C) , (46) 



and Ai,2(M) and Si,2,3(k) should be taken equal to zero at k such that ei > 1, i.e., when the following inequality 
satisfies: 

1 - /3 - 3C < 0. (47) 

Notice that £(g,k) = at g = ±^/l~^~ei, \/T~\~e2 and it is positive inside the intervals (— -\/l~^77, vT^^^e^) and 
{VTT7^,oo). When P,^ < 1 one has Ai{k) w - ln(/3 + VC), -42.3(k) - l/(/3 + V?) and Si,2,3(k) « 1. If (3,£,>1 
we have .43(k) ~ l/V/3 + Terms in Eq. (|36p containing ;Bi(k) and B4 play only at | sin2(/)k| <C 1. In particular, we 
obtain from Eq. ([36]) that the spin- wave damping is zero at T = and such k that £1 > 1, i.e., when inequality (|47p 
holds. The region determined by Eq. (|47|) is sketched in Fig. H) 

We conclude from Eqs. and ([55)1 that the damping is small in the corresponding intervals provided that 



C. k-^Scjo/D 

As it is found above, the real part of the spectrum at fc < {A/(Sujq))'^ is renormalized greatly being equal approxi- 
mately to A. We find self-consistently for Suq/D > fc > A/VtD 

^■^^2^-05^- 

Notice that this regime is realized at large enough temperature, T '3> ujq. At smaller k, k <^ A/^/TD <^ Sloq/D, or 
at small temperature, T <ti luq, the spin- wave damping is exponentially small: 

r.ocexp(-^). (49) 

It should be noted from Eqs. ([55]) and ([^5)1 that when T 3> 6*1^0 the damping Ft increases upon decreasing k for 
k < {Sujo/Df/^. On the other hand, as we have just obtained, the damping is exponentially small at very small 
momenta k <C A/^/TD. Hence, one should observe a peak in Fk at A: /S./\/TD and at any k direction which height 
can be estimated from Eq. ([48]) . In particular, if the temperature is as large as the interval {/S./ {Su}q)Y 3> fc > A/ \/TD 
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is finite, i.e., if T ^ S'^ujq/C, we fiave ek ~ A at fc ~ A/\/TD in Eq. ([48]) and one obtains for the peak height from 
Eq. mi) using Eq. ^ 

In general case the peak height cannot be larger than the value given by Eq. ([50)) because et > A at fc ^ Sujq/D. 
Notice also that Eq. ([50)) is valid for arbitrary (/>k. 

It is seen from Eq. (|50p that the spin- wave damping is much smaller than the real part of the spectrum if T ^ Z?, Tc, 
where 

ln{47:S{D/[SLUo]r/^) ^ ' 

is the value of the Curie temperature for 5 ~ 1 obtained using the spin- wave theory (see Sec. l Vip i^ On the other hand 
the temperature can be greater than D for large spin values much greater than unity, 5 > In (47rS'(Z?/[S'(jJo])'^^^) 
so that D <T<^ T^^^'^\ where 

(s»i) _ 8nDS 

- 3HD/[Su;o]) ^''^ 

is the Curie temperature in the spin-wave approximation for S S> In (47rS'(£'/[S'(jJo])'^^^) ■ Then one might conclude 
from Eq. ()50[) that the damping can be much larger than the real part of the spectrum. Meantime we show in the 
next section that the temperature correction to the spin-wave gap is large at such large T and S. As a result the 
imaginary part of the spectrum is also much smaller than the real part at S" 3> 1 and D < T <^ j.OS»i) ^-^^ 
peak height of the ratio Fk/ek m quantum 2D FM cannot be larger than that in the classical 2D FM which is equal 
approximately to 0.16 for the simple square lattice with exchange coupling between nearest spins only. 

We sketch the dependence of Fk on the momentum in Fig. |5|atfc^l,5^1 and T ^ SuJa taking into account 
the results obtained in this section. It is seen that the damping is highly anisotropic at fc > Sluq/D. The damping 
increases with decreasing k up to k ^ A/ \/TD if k is directed along a square edge (i.e., if | sin 20k | =0). In contrast 
the damping is not monotonic function of k for | sin20k| ~ 1: it decreases with decreasing fc up to fc ^ (SuJo/Dy^'^ 
and then it rises up to fc ^ A/\/TD. The damping is only slightly anisotropic in the interval A/^/TD ^ k < Sujq/ D. 
There is the peak aX k A/ \/TD at any given (/)k which height can be estimated using Eq. ([48)) . This peak is followed 
by exponential decay of the damping at fc < A/^/TD having the form \4Q\ . 

To illustrate in more detail the region of small momenta k ^ Scuq/D in which the peak exists, we derive general 
expressions for Fk which coincide with Eq. ([36]) at /c > Sujq/D and with Eqs. ([48]) and ([50]) at A/VtD < fc <C Sloq/D. 
These general expressions appear to be quite cumbersome for arbitrary 0k- Then, we present here only the equation 
in the special case of sin 0k = which is the most simple one: 

exp(gVl + <;Vi) g(i+g^f/^ /, l + q' 1 + k' 
''(e.p(,yTT7/<)-iy<l + WV + ' '''' 



Fk _ _c^S\j^iVTTk_ 



C 



where k = ky/SDujQa/ A, t = T/[SbjQa) and 

C - (l - 3k2 + V9«;4 + 10^2 + l) . (54) 

At q — C, the expression under the square root in Eq. ()53p is equal to zero and it is positive for q > C,. When 
1 > K > ^/SuJoa/T (i.e., when A/VS'Dwo > fc > A/VtD) and T > Sujo, Eq. ([53]) transforms into Eq. HH)). In the 
opposite limiting case of k ^ 1 (to be precise, at A/ ^ SDujq <C fc ^ -y/ Sluq/D) one obtains the last term in Eq. ([36)) 
from Eq. ([55| at T » ^^n. 

We plot in Fig. [6] the ratio of the spin- wave damping and the real part of the spectrum given by Eq. (|53p versus the 
reduced wave-vector k for 2D FM with 5 = 1/2 and S* = 3 on the simple square lattice assuming that ujq = O.OIJ. 
The peak is seen at k ^ SuJoa/T (i.e., at fc ~ A/\/TD). Its position moves to smaller k and the height rises as S 
increases at a given ratio T /Tc or as T increases at a given S. 
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V. LARGE SPINS 



Let us consider now large spins 5 > In (47r5(-D/[5ti'o])'^^^) • As it is explained in the previous section, the temperature 
can be of the order of D in this case remaining much smaller than the Curie temperature. When T D one can use 
expressions for U,{u), k) obtained above with the exception for the gap: the temperature correction becomes important 
and one has for the gap in Eq. ([32]) 



SujQa 



q 



1 + ^ (1 + iVq) 



(55) 



where A is given by Eq. (|3ip . the first and the second terms in the square brackets stem from the Hartree-Fock and 
the loop diagrams shown in Fig. [2lb) and (c) , respectively. 

Let us discuss temperatures 25 Jo <T <^ T^^^\ where 2SJq is the spin- wave band width and T^^^^ is given by 
Eq. (|52p . To find corrections to the spectrum at such large T one can expand all Plank's functions under sums over 
momenta up to the first term: TVq « T/e^. Then, we have Eq. ([5^ for the real part of 0(tt>, k), where 



V{T) 



2ttD "^KSujo 



(56) 



now and one obtains for the gap from Eq. ([55l 



A» = ^/A2+aC»52-^T, 



C» = 



2f2 



(57) 
(58) 



where we can discard A^ under the square root in Eq. ([F7|) because of its smallness compared to the second term at 
T ^ D. Summation over large momenta gives the main contribution in Eq. (|58p and we find as a result of numerical 
computation for the coupling between only nearest neighbor spins on the simple square lattice C> ~ 0.025. 

A new wide region appears in the momentum space in which the bare spectrum is much larger than the gap and has 
the form el « SuJoa{DP + {SuJok/2) sin^ c/)^ ). The corresponding interval in the k-space is given by (A^ / (Sluq))'^ < 
k ^ Sivo/D for I sin 0k I ^ 1 and by A> / V SDujq ^ fc ^ Sluq/D for |sin(/)k| ^ 1. As a result all the expressions 
for the damping obtained above should be reconsidered. In particular, Eq. is valid only for k ^ Sy^ujo/T. All 
T-independent terms are negligible in Eq. (|34p and one should use expression (|57p for the gap rather than Eq. (|3ip . 
Then, 



Eqs. (1211) and are not valid at ^JSoJo/D < fc < S^Juoo/T and Sujo/D < fc < ^JSuo/D, respectively, 
because integration over the above mentioned interval in the momentum space is essential in the corresponding sums. 
The results are quite cumbersome and we do not present them here. The main conclusion is that the spin-wave 
damping is much smaller than the real part of the spectrum at 3> Sujq/D. Let us discuss in somewhat detail only 
the damping at small momenta k <C Suiq/D because the diffusion mode was proposed in Refs.SiiS at small k. 

The ratio of the spin-wave damping and the real part of the spectrum at {A-^ / (Sujo))'^ ^ k > A-^/^/TD can be 
found from Eq. ((50)) (replacing A with A^,) with the result 



A 



1 



> 



28 a 



(59) 



> 



This expression gives the approximate value of the peak height on the curve Fk/ek that is valid for all t/ik. The 
right part of Eq. (|59p is equal approximately to 0.16 for the exchange coupling between only nearest neighbor spins 
on the simple square lattice. It is interesting to note that the smallness of the spin-wave damping in this case is 
numerical whereas at 5* ~ 1 the smallness is parametric (see Eq. ([SO)) ). Notice that at fc ^ A^/^/TD the damping 
is exponentially small, as it is discussed in Sec. IIV CI 

It is seen from Eq. ([55]) that the value of the peak height is in inverse proportion to which depends on the 
exchange coupling between spins Jj^ (see Eq. (jSSp ). It is interesting to examine the dependence of Fk/A^ given by 
Eq. (|59p on the value of the coupling between next-nearest neighbors which, in particular, can reduce D significantly. 
We assume that the exchange coupling between nearest- and next-nearest neighbor spins are equal to J and J', 
respectively (see the inset in Fig. [7]). One has for the spin-wave stiffness in this case D = SJ{1 + 2J'/J). The 
dependence of the peak height on D/{SJ) is shown in Fig. [71 It is seen that even at very small D, i.e., for frustrating 
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next-nearest-neighbor interaction, magnons are well-defined quasi-particles. Importantly, it is implied even in the case 
of smah D that D » Suq and 8S(J J') <C T < SD, where 8S{J + J') is the spin-wave band width. 

We do not present here expressions for the damping at {A^ / {Sluq))'^ ^ k < Suo/D for arbitrary (pu. The main 
conclusion is that the damping is much smaller than the real part of the spectrum. To illustrate this let us consider 
only sin(/)k = 0. The bare spectrum has the minimum value at sin(/)k = for a given k in the discussed momentum 
interval and it is the most "dangerous" case in which one could expect large spin- wave damping as compared with 
the real part of the spectrum. Corresponding calculations give Eq. (|53p for the damping at fc ^ Sujq/ D, where one 
should replace A with A^,. We use this modification of Eq. (|53|) to plot Fk/ek in Fig. [6] for S = 30, ojq = O.OIJ and 



Classical spins 

To find the spectrum renormalization in the classical 2D FM using the spin- wave formalism discussed in the present 
paper one should consider the limit of 

5^00, h^O, J,wo^O (60) 

assuming that 

hS = const, JS^ = j = const, ujqS^ = w = const (61) 

and T/j is much smaller than unity. Moreover, one should replace operators of creation and annihilation ak, with 
classical operators /3k = Ok/v^, /3k = al^/VS which Bose occupation numbers are finite.— As a result the spectrum 
and corrections to it in the classical limit can be obtained from the expressions found above by multiplying them by 
S and taking the limit ((60|) with the assumptions ([61]) . In particular, we have for the gap in the classical 2D FM from 
Eq. jSZl) 



A^ = ^aCy>-jT. (62) 

Only the second term under the square root in Eq. ([57]) contributes to A^ . Notice that there is no exponential decay 
of the damping at small k in classical 2D FM in which Eq. ([59]) is valid at all k much smaller than C^wT/ p because 
A^/VtD — > when S oo. Then, the peak in Fk and Fk/ck is located at fc = 0. 

In particular, one obtains in the limiting case of classical spins from Eq. (|53p by expanding the exponents in t and 
replacing A with Aqo 



Fk 1 VTT^ r , (1 + 9')'/' L 1 + 9' 1 + ^:' ,a^^ 



ek 167rC» k3 J (■ ^g(l-h2g2)3y (1 + 2^2)2 ^2 



where k, = k^J p / [C^w'^T) now and Q is given by Eq. ([51)) . Eq. transforms into Eq. (|59p at k ^ 1 (with Aqo put 
instead of A^). We plot in Fig. [6] also the ratio of the spin- wave damping and the real part of the spectrum given by 
Eq. ([63]) versus the reduced wave-vector n for the classical 2D FM on the simple square lattice. As it is pointed out 
above, the position of the peak on quantum curves moves to smaller k and the peak height rises as S increases at a 
given ratio T /Tc or as T increases at a given S. But as it is demonstrated above and as it is seen by the example 
of 5 = 1/2, 5 = 3 and 5 = 30 shown in Fig. [6l the peaks on the quantum curves cannot be higher than that of the 
classical 2D FM located at fc = which height is given by Eq. (f59|) . It should be noted here that the peak height 
for finite 5^1 and D ^ T <C T^'^^^ is slightly smaller than that given by Eq. ([55)1 because we have discarded A^ 
under the square root in Eq. (|57p deriving Eq. ([5^ . On the other hand Eq. (|59p gives the precise value of the peak 
height for classical 2D FM because A^ disappears after taking the limit (|60|) with the assumptions (|6T|l and only the 
second term under the square root in Eq. (|57p contributes to the gap. 



VI. DISCUSSION 



We address in this section five questions: (i) relation between the anisotropic term ([^0)1 in the total energy of 2D FM 
and the spin- wave gap pip , (ii) discussion of the further order 1/5 corrections, (iii) calculation of the magnetization 
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taking into account the spin-wave spectrum renormalization, (iv) derivation of the spin Green's functions in the first 
order of 1/5*, and (v) brief discussion of the effect of the easy-plane anisotropy and the spectrum renormahzation 
obtained in Refs.— i^^. 

(i) As is discussed in Introduction, it looks reasonable that the dipolar in-plain anisotropy given by Eq. (j20p should 
be accompanied with a spin-wave gap. The relation between this anisotropy and the gap can be shown for 2D FM 
in the same non-rigorous manner as for 3D FM.— Let us discuss large spins and try to take into account the dipolar 
in-plain anisotropy (|20p phenomenologically by adding to the microscopic Hamiltonian ^ the following expression 
[cf. Eq. (HOI)]: 



^^E(^r)'(^f)'- (64) 

This term after Dyson-Maleev transformation ([5]) gives the contribution CSujQ/{4:D)(al^ + Ck)'^ to the bilinear part 
(HI) of the Hamiltonian that in turn leads to the shift E'k >-> i^k + C Soul/ {2D) and B]^ i-^ + CSujI/{2D). Using 
this renormalization of ii'k and i?k and Eq. for the spectrum one recovers the spin- wave gap (pij) in Eq. (|30p . 
This consideration does not work for S = 1/2 because Eq. (fM]) is a constant in this case. 

(ii) Let us turn to further order 1/5 corrections. Some diagrams of the second order in 1/5 are presented in 
Fig. [9l It can be shown that at least their real parts are much smaller than the real parts of the first order diagrams 
discussed above. To demonstrate this one has to take into account that the spin-wave gap screens infrared singularities 
appearing in some of these diagrams. Moreover three- and four-particle vertexes contain additional smallnesses at 
small external momenta. Thus, the three-particle vertex (|10p contains a;2;-component of the dipolar tensor which is 
proportional to the product of loq and momentum (see Eq. ([5])). As for four-particles vertex, the expression under the 
sum in Eq. pTjl has the form at fci.2,3.4 ^ 1 

2k3(2ki -I- k3)alial2a3a4 - 
l + sin2(/)k2 t , 

^ a„;^a2a3a4 + (65) 

^ (k2 sin^ + 2 |k2 + k3| cos^ ^ka+k:,) a^ia!_2a3a4 + 
aujooLi (a-z — ) 0304- 



It is seen that the first term in Eq. (j65p is quadratic in momenta. The second and the third ones are proportional to 
and momenta. The last term in Eq. (|65p is proportional only to ujq but the combination 02 — 0^2 involved in it is "soft" : 
for instance, it's coupling with operator a_2 gives F{uj2,'k2) — G{cj2 7^2) ~ — (DA:| -f 5wofc2 sin^ 0k2/2 — a;2)/2?(<^2, k2) 
that is much smaller than F{uj2,'k2) or G(a;2,k2) themselves (F(w2,k2) « G'(w2,k2) ~ Sujo/'D{uj2,^2) at k2 <C 
yJSiJi)/!) and uj2 ^ 5aJo)- As a result further order diagrams appear to be small at T ^ Tc- 

Unfortunately, the smallness of the three- and four-particle vertexes was not taken into account in my previous 
paper— devoted to 3D FM in similar qualitative discussion of further order 1/5 corrections. As a result further order 
corrections were overestimated there. It was proposed that they are small only at T ^ wq, whereas the range of the 
validity of the perturbation theory is much wider in 3D FM: T <C Tq. 

(iii) Let us calculate now the magnetization value (5^) using renormalized spectrum. We have after simple compu- 
tation 



5 91 ^ 25eq I ujoa Sujq 

I 167rL>' ln(i:>/[5wo])' 

Notice that zero-point fluctuations give the main contribution at T ^ 5a;o/ lii(-D/ [5^^)01 )■ The value 5 — (5^) given 
by Eq. ([66]) is shown in Fig. [8] for luq = O.IJ. It should be pointed out also that taking into consideration the gap 
in the spectrum does not change the form of the magnetization. Meantime Eq. (|66p differs from the corresponding 
expression in Ref .-. The origin of this discrepancy is discussed in Appendix |B] The value of the Curie temperature 
([5T|l giving by the spin-wave theory follows from Eq. ([TO]) by putting (5^) to be equal to zero. In the case of large 
spins 5 > In {4:TTS{D/[Sujo]f^^) the Plank's funct ion in Eq. (|66p can be expanded and one finds Eq. (j52p for the Curie 
temperature by putting (5^) = 0. 

(iv) Spin Green's functions defined as Xij(w,k) = i dte^'^'^ {[S^_j^{t) , Sl_{0)]) can be calculated straightforwardly 
in the first order of 1/5 using Eqs. (|16p - (fl^ and expressions for the self-energy parts with the following results for 
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the transverse components 



, ^ Dk^ + (^c^o/2)fcsin^ cf,^ + AV(5a;oa) 



where Ck is the bare spectrum here. Corrections from the self-energy parts to numerator are negUgible in Eq. (|67p . 
In contrast numerator in Eq. (j68p for Xyyi'-^i^) renormaUzes greatly at small k so that the uniform susceptibility 
Xyyi(^,0) becomes finite. 

The corresponding expression for the longitudinal component Xzzi^^,^) is slightly more cumbersome. It can be 
calculated using Eqs. (|A2p . Let us discuss only uniform longitudinal susceptibility Xzz{^j->,0)- We have for ^ 
^{Sioor/D » A 

2A 



/ n^ TSujoa 1/2 
One obtains for the imaginary part when \uj\ < a/ {Su!o)^/D 



isgn(Lj) I . (69) 



lmx..(c.,0) =sgnH0(|^| -2A) ^^3/4 ^2 ^ (^0) 



where r(a:) is the gamma-function and 6{x) is the thcta-function (^?(a;) = 1 when a; > and 0{x) ~ when a; < 0). 
The corresponding expression for the real part of Xzz{'^^^") is quit cumbersome and we do not present it here. It is 
seen from Eq. (|70p that Imxzz(a;, 0) = at w < 2A. In contrast, if one did not take into account the spin- wave gap 
A the infrared singularity would appear of the form Imx2z(a; — > 0, 0) ^ lo^'^I'^ . Such a singularity is nonphysical one 
since it leads to an infinitely large absorption function^ Qt^ oc ti;Imx(cj) at = 0. Similar situation exists in 3D FM 
with dipolar forces. The infrared singularity of the form \varxzz{}^ — > 0, 0) ~ T ju can be found not taking into account 
the spin- wave gapi^ This singularity leads to the finite absorption function at w = signifying that the sample would 
be heated by a dc field. On the other hand the spin- wave gap screens this singularity leading to zeroth Q^^ at w = Ofii 
(v) We will discuss in detail elsewhere the effect of the exchange anisotropy having the form 



na = \Y^M^s\si,. (71) 



and the one-ion anisotropy The corresponding rcnormalization of the above results for the real part of 

the spectrum is discussed briefly in Appendix [C] taking into account the exchange anisotropy (171]) . In particular, it is 
shown there that expressions (|3T|) and ((57| for the gap should be multiplied by the factor ^/Gm^Juoq, where 

wo=t^o + — • (72) 
a 

In the limiting case of classical spins one should imply in addition to Eq. (|60p and Eq. ([6T|) that ^4 ^ and AS^ — const. 
Then, we will assume that uj^S"^ — w — const and one has for the gap in classical 2D FM Eq. ([62|) multiplied by 



Classical 2D FM with dipolar interaction and easy-plane anisotropy was discussed in Refs.^'^". Exi sten ce of the 
easy-plane anisotropy seems to be not crucial for the results obtained in Refsi^ii^. The condition \/w/ j 3> w/j, 
that is important for the consideration in these papers, holds also aX A = Q (w = w) \i w <^ j that is also implied 
there. The spin-wave gap ([S^ (multiplied by ^w/w) was not taken into account in Refsi^i^^. Great renormalization 
of the spin-wave spectrum at small enough momenta was obtained there. In particular, a diffusion mode was found 
at fc <C koM ~ wt'^/'^ /\j \v}/^{yj jw /w)], where t = T/{Anj). Meantime the spin-wave gap screens all the spectrum 
peculiarities obtained in Refsi^iiS. For example, at jsin^kl ^ 1 the energy of the diffusion mode has the form up 
to a numerical factor of the order of unity — ifc^i"^/"* In^^^ {\/ jw /w) ^ /w"^ . The spin-wave gap given by Eq. 
multiplied by ^Jw/w is much larger than the energy of the diffusion mode at fc ^ k^M- 



VII. SUMMARY AND CONCLUSION 



In the present paper we discuss two-dimensional Heisenberg ferromagnet (2D FM) with dipolar forces at < T <C Tc 
described by the Hamiltonian ^ and consider renormalization of the bare spin- wave spectrum p4p due to interaction 
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between spin waves. For this purpose we carry out a comprehensive analysis of the first 1/S' corrections to the 
spin-wave spectrum originating from diagrams shown in Fig. [21 

We obtain the following results for S" ~ 1. Corrections to the square of the real part of the spectrum are given 
by Eqs. ([50]) and ([5^ . In particular, at T ^ Sujq, where ujq is the characteristic dipolar energy given by Eq. ([1]), 
thermal corrections result simply in renormalization (j28p of the constants D (spin-wave stiffness) and luq in the bare 
spectrum ([T4l) . This renormalization is small at T ^ Tc- But it is significant that similar to 3D FM considered in 
our previous paper^-'^ we obtain also the spin- wave gap A in the spectrum given by Eq. (|3ip . This gap stemming from 
the spin-wave interaction renormalizes greatly the bare gapless spectrum (ITil) at fc < {A/ (Suiq))'^ (see Fig. [3]). 

Spin- wave damping Fk at T = is given by Eqs. ([M]) and ([55)1 for momenta 1^ yj Sujq / D and k <C \/SuJo/D, 
respectively. There is a region in the k space at small momenta sketched in Fig. [31 in which the damping is equal to 
zero. Meantime we find great thermal enhancement of the damping in this region. 

Temperature fluctuations give the main contribution to the spin-wave damping at T 3> Sloq that is given by 
Eqs. ([SI, dSi), ([48]), and ^ for momenta k » ^JSloq/D, Sluq/D <C fc <C y/SuJo/D, A/VTD < fc < Sluq/D, and 
fc < A/VTD, respectively. The dependence of the damping on the momentum is sketched in Fig. [H for fc <C 1. It is 
seen that Fk is highly anisotropic at fc > Sujq/D. The damping increases with decreasing fc up to fc ^ Aj \JTD if k 
is directed along a square edge (i.e., if | sin20i(;| =0). In contrast the damping is not monotonic function of fc when 
I sin20i(;| ^ 1: it decreases with decreasing fc up to fc ^ [SloqI D)'^!'^ and then it rises up to fc ~ A/^/TD. The damping 
is only slightly anisotropic in the interval A/\/TD < fc < Sloq/D. There is a peak at fc ~ A/\/TD at any given 0k 
which height can be estimated using Eq. (I48p . If the temperature is large enough so that the interval is finite given 
by {A/ {SuJo)Y ^ k > A/^/TD, the peak height is given by Eq. ([50)1 . This peak is followed by exponential decay of 
the damping at fc < A/ \JTD having the form ([^5)1 . 

An important quantity to be examined is the ratio of the spin-wave damping and the real part of the spectrum 
Fk/sk which is much smaller than unity if magnons are well-defined quasi-particles. It is seen from Eqs. p4p . (j36p . 
and (|48p that the ratio Fk/ck rises upon decreasing fc at fc > Sujq/D for all 0k- This growth changes into exponential 
decay at fc < A/^/TD given by Eq. ([49)) . Then, there is a peak in Fk/ck at fc ~ A/^/TD and at any given 0k- This 
peak rises and its position moves to smaller fc as S" increases at a given ratio T /Tq or as T increases at a given S (see 
Fig. [6]). Meantime its height is restricted by Eq. ([50]) which is proportional to T/D. In the case of 5* '--^ 1 we have 
T/D < 1 when T < Tc- 

On the other hand, as it is pointed out in Sec. IIVI C, the temperature can be greater than D for large spin 
values much greater than unity, S* > In {4,-kS{D /[Sloq])^/'^^ , so that D < T <^ Tc (x DS. Renormalization of the 
spectrum should be reconsidered at such large T and S because new large temperature 1/ S corrections arise. Such 
reconsideration is carried out in Sec. [Vl In particular, we find thermal renormalization (I28p of the constants D and wq 
in the bare spectrum ([T4l) . where W{T) and V{T) are given by Eqs. (I56p in this case. Then, thermal correction to the 
gap becomes important and we have Eq. ([55)1 for the gap sAT ^ D which transform into Eq. ([57)1 at 25* Jo ^ T <C lb, 
where 2SJq is the spin-wave band width. Focusing on the spin-wave damping renormalization at small fc only, we 
observe that the peak height of Fk/ek cannot exceed the value given by Eq. (|59p that is equal approximately to 0.16 
for the simple square lattice and that is a counterpart of Eq. (jSOp for 5^1. It is interesting to note the numerical 
smallness of the peak height in Fk/ck when S* ^ 1 and Tc ^ T ^ 2SJo that contrasts to the case of small 5*, in 
which the peak in Fk/ck cannot exceed the value ([50]) proportional to T/D <C 1. The small value given by Eq. ([59]) 
increases only slightly upon taking into account a frustrating ncxt-nearest-neighbor exchange coupling (see Fig. [7]). 
The limiting case of classical spins is also discussed in Sec. [Vl at T ^ j, where j is the exchange constant in the 
classical model (see Eq. ([61]) ). We obtain expression ([62]) for the gap and the peak in Fk/ek at fc = which height is 
given by Eq. (|59p precisely (see Fig. [6]). Thus, we find that magnons are well-defined quasi-particles in both quantum 
and classical 2D FMs with dipolar forces. 

We note that appearance of the gap is accompanied by the anisotropy in the total energy of the quantum 2D 
FM given by Eq. ([20]) and caused by quantum fluctuations lifting the degeneracy of the classical ground state. We 
demonstrate in Sec. IVII the relation between the dipolar anisotropic term (|20p in the total energy and the gap in the 
spectrum at T = 0. 

Spin Green's functions Xiji'^i'^) are derived in Sec. IVII with the results ([57]) and ([55]) for the transverse diagonal 
components Xxx{^i k) and Xyy{^^ k), and (|69p and ()70p for uniform longitudinal one Xzz('^, 0). It should be noted that 
if one did not take into account the spin- wave gap the infrared singularity would appear of the form 1tiixzz{^ ^ 0, 0) '--^ 
uj^'^/'^. Such a singularity is nonphysical one since it leads to an infinitely large absorption function Q^^ oc ojlmxito) 
at cj = 0. The spin- wave gap screens this singularity: as it is seen from Eq. ([70]) . Imxzz(ti',0) = at w < 2A. 

Modification of the results by taking into consideration the easy-plane exchange anisotropy ([7T]) is discussed briefly 
in Sec. IVII In particular, we have shown that expressions ([3T]) and ([57]) for the gap should be multiplied by the factor 
i/ifjoTwo, where wq is given by Eq. ([72]) . 

Expression (|66p for the magnetization in 2D FM is obtained which differs from the well-known result of Ref.— . 



15 



Higher order corrections to the spectrum are discussed in Sec. I VII and it is concluded that they are small compared 
to the first corrections obtained. 

We would like to note in conclusion that the spectrum is gapless in 3D antiferromagnets with dipolar forces in the 
spin- wave approximation^ and the spin-wave interaction should lead to the gap in 3D antiferromagnets similar to 2D 
and 3D FMs. 
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APPENDIX A: CALCULATION OF n^^\uj,k) 

We present in this appendix some details of calculation of Q^"^") (uj,'k) that is a contribution to il(aj,k) given by 
Eq. (I18|) from the loop diagram shown in Fig. [2jc). This diagram originates from Ti.^ terms (|10p in the Hamiltonian. 
As a result of simple calculations we lead to quit a cumbersome expression 



xi^iQ^^' f {B1B2 + E1E2+UJ1UJ2) (Ala) 

+ 2Qi^^Qr([Si - B,][E2 - B2] + c^ic^2) (Alb) 

+ Ql'QT{[Ei - Bi][E2 - B2] - c^it^2) + 2 {Ql' f Ei{E2 - ^2)) (Ale) 

^ ^f-^_^[{^L0,f-el][{^U2)^^el] 

X (qY{QT + QT){Ei - B,){E2 - B2) (Aid) 

+ [Ql'f {B1B2 + E1E2 + UJ1UJ2) (Ale) 

+ 2Ql'Qr{[Ei - B,][E2 - B2] + iJiuj2)) (Alf) 
25 X ^ iui 



[q1^qT + QT){E2-B2) (Aig) 

- Ql'[QT{E2 + B2) + or (^2 - E2)]) , (Alh) 



where k = (k, cj), gi 2 = (qi.2j'^i.2) and we drop the index q in Eq. (jAip to light the notation. Sums S + E and 
S + E + n + n^ lead to terms ((XTb )- ((XTb ) and ((XIH )- ((Xll "). respectively, whereas terms (RTk ) and (PIl i) result from 
iw(E-E) (see Eg. ([TSl)). 

One obtains the following expressions after summation over imaginary frequencies and analytical continuation on 
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Lo from imaginary axis to the real one: 

1 



l + 27V(ei) , l + 2iV(e2) , 2{ei ~ e2){N{e2) - N{€i)) 



2ei[(ei+e2)2-(u; + z5)2] 2e2[(ei + £2)^ - (c. + *J)2] [(e, + e^^ - {cu + zSm^, - 62)' - [to + iS)^ 
{iLUi){iLU — iuji) 



■(A2a) 



ei(l + 27V(ei)) , e2(l + 27V(e2)) 2e,e2iei ^ 62^(62) ~ N{e,)) 



(A2b) 



2[(ei + £2)2 -{u; + tS)^] 2[(ei + £2)2 - {uj + tS)^] [(ei + £2)2 - {tu + tS)^][{ei - €2^ - {uj + i5Y] ' 

( l + 2N{e2) 2ei(Ar(e2)-jV(6i)) \ 

\2e2[{ei+e2Y -[uj + iSf] ^ [{e^ + e2Y - + i5Y][{ei - e2f - + i5f] ) ' 

We calculate now the real part of ri^'^)(aj, k) using Eqs. (|A2|1 . (jAld ) is the only term remaining finite at aj,k = 
and leading to the contribution ([29]) to the spin- wave gap and to the second term in the square brackets in Eq. (j57p . 
The rest corrections in Eq. (|A1[) are much smaller than either (jAld ) or corrections from the Hartree-Fock diagram 

given by Eqs. ^ and (US]). Let us estimate them. Term (|XTb ) is of the order of ui^^^k^^^ sin^ 2(j)k{ky^uJaD + T) /^/D, 
wgfc^sin^ 2(j3]^{ky/SujQ/D + T / D) and cjg sin^ 2(/)k(wg + k^TD)/D^ at /c < Sujo/D, Sujq/D < A: < ^SuJn/D and 
k ^ ^ Suiq/D, respectively. Term (|Alb ) is of the order of fc2ci;g. It is much smaller than the first term in Eq. ([24|l if 
ln(Z)/<x'o) ^ 1 that we assume to be held. Contributions (|Alb ) cannot give more than (Klh) and (|XTb V Term ((JTb ) 
is of the order oi uj^^'^Vksm'^ 2(/)k(fc\/^^ + ^)/-D^/^ wgsin^ 2(t)u{k^/U^ + T) / D'^ and wj^ sin^ 2(j)y,{u}l / k"^ + T D) / D'^ 
at fc < Sll>o/D, Sujq/D < fc < y^SuJo/D and fc > ^JSujq/D, respectively. ([JIT :') does not exceed (IXIb V (IXTk ) is of 
the order of Ld"^ [luI / D'^){\ii{D / ujq) + T/(wo + Dk"^)). jJlJi) does not exceed the sum of (|XTb ) and (|XTf). 

Let us turn now to the imaginary part of US'^\uj, k). Corresponding calculations have been done straightforwardly 
using Eqs. (jAl[) and (|A2[) with the following results. 

fc > ^Sujq/D. Terms (|XTt) and (PIk ) give the main equal contributions at small temperature leading to T- 
independent terms in Eq. ([34]) . Terms (jAlb ) and (|Al^ ) give equal contributions to the last term in Eq. ((34|) . The 
second term in the square brackets in Eq. ([34]) containing /(0k) originates from (|Alh ). (lAlb ) and (jAlb ): endowments 
of (jAlb ) and (jAlb ) are equal and twice as little as that of (jAlb ). The origin of the function /(0k) is the following. Sum 
of the form 0T~^ (5(cos(0k — 0q) — eq/ C^Dkq)) / {qe'^) appears in the corresponding expressions in which summation 
over small momenta q <C Sloq/D is essential. We have for the square of the spectrum at such q (cf. Eq. (j55|) ) 

e^«^^(.2+.sm2 0q + ^j, (A3) 

where q = q2D/ (Sloq). 

Sojq/D <^ k <^ ^ Sloq/D. One finds Eq. (|36|) in this regime. Term in Eq. (j36p containing -4i^2,3(q) comes from 
(jAlti )~ (jAlf ) whereas those containing y8i(k), S2(k) and By,{\<i) originate from (jAlb ). (jAl^ ) and (jAld ). respectively. 
Term in Eq. (j36p proportional to B4, comes from (lAlt i)- (IAlF ). 

We demonstrate now in what way the quantity £(g, k) appears in Eq. (j34|) . When Sujq/D <C fc ^ ^/SLdo/D it is not 
sufficient to use the leading term in the expression for the spectrum assuming that ek ~ \fSaujoDk. Really, w e have in 
this case (5(ek-eq-ek+q) oc ^((^k-'/'q-Tr)^) and J(ek-eq+ek+q) oc ^((^k-^q)^) if Sluq/D < k, q, |k + q| < ^JSujq/D. 
At the same time expressions under the sums do not vanish at (j}^ — or — n. Then, the appearance of the 
squares in arguments of delta-functions signifies that one should take into account smaller terms in the expression for 
the spectrum: 

.,.V^(,.^^s.„'^, + ^ + ^). (A4) 
Using Eq. (|M| we have if Sluq/D < k, q, |k + q| < ^JSiOalD 



ek - e, ± ek+q « ^^(k^) ['^^ - 16^^ ('a: - ^'^) ) ' 
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where q_L(q|| ) is the component of q perpendicular (parallel) to k. We assume in Eq. (|A5p that — jgy | ± |fc + | = 
and q± <C , k. 

k ^ Sujq/D. In this regime Imfi'^' (w, k) is finite at high temperature only and one has Eq. (j48| for it. Term ()A1H ) 
only contributes to Eq. 



APPENDIX B: DISCUSSION OF THE DISCREPANCY BETWEEN EQUATION ([MD FOR THE 
MAGNETIZATION AND THE PREVIOUS RESULT 

In this appendix we comment on the discrepancy between Eq. (ITOj) for the magnetization and the corresponding 
expression in Ref."*. Eq. coincides with that obtained in Rcf."* at T ^ Sloq /ln{D / [Scoo]) and at T = 0. Meantime, 
a term was obtained in Refi^ proportional to T^^^ujq at T ^ cjq^^/Vj that is much larger than the T-independent 
term in Eq. (|66|) at loq^^ /VJ ^ T ^ J{luq/ J)^^/^ . I believe that the term proportional to T^/^Wq observed in Refj^ 
is an artifact. Careful calculation with the bare spectrum ([Ti|) gives the term proportional to T^/^uj^^^^ rather than 
T^^'^luq at T ^ Wq^^/\/J which is much smaller than the T-independent term. To confirm this finding we note 
that the values of (S — {S^-))/S obtained in the regimes of T > uil^^/Vl and T < luI^^/v^ should be of the same 
order at T ~ ojq^^/a/J. At the same time T^^'^ujq '^^^ ^ ^/uJ^ and T\vl{T juf'^'^) ^ u?^"^ at T ^ ujf/^ Surprisingly, 
this artifact has not been revealed so far (see, e.g., Refsi^!^). To conclude, one leads to the same expression (|66p for 
{S — {S^))/S as a result of calculations using the bare (fH)) and the renormalized spectra. 



APPENDIX C: EFFECT OF THE EXCHANGE ANISOTROPY 



We discuss briefly in this appendix the effect of the exchange anisotropy given by Eq. (|7T|) . We consider here only 
the exchange anisotropy which differs from the on-site one having the form A'^i{Sf)'^ that exists in thin ferromagnetic 
films. 1 The reason is that it is technically easier to discuss the exchange anisotropy (|7T|) . Moreover it is believed that 
these two types of anisotropies lead to similar physical results (see Refi^ and references therein). A more detailed 
discussion of the effect of the anisotropy in 2D FM with dipolar forces will be published elsewhere. We consider below 
both signs of A;,„, i.e. both easy-axis and easy-plane anisotropy. Easy-axis one competes with the dipolar anisotropy 
which favors in-plane spins alignment. We restrict ourself here to the case of not too large easy-axis anisotropy at 
which spins lie within the plane. 

Adding Ha to the Hamiltonian ^ one obtains after the Dyson-Maleev transformation ([5]) a renormalization of i?k 
and i?k in the bilinear part of the Hamiltonian ^ 



(CI) 



and a contribution to the four-magnon term 



n 



(a) 



1 

4^, 



(C2) 



ki+k2+k3+k4=0 

As a result of renormalization (|C1[) the bare spectrum has the form at A; <C 1 (cf. Eq. HH)) 



' [Dk'^ + SCo^a) ( Z?fc2 ^ sin^ , 



(C3) 



where Qq is given by Eq. ([72]) . Notice that in the case of easy-axis anisotropy {A < 0) the spectrum e^^"^ becomes 
imaginary at small enough k if ujq < 0- At the same time the in-plane spin alignment becomes energetically unfavorable 
if the anisotropy is as large as (Dq < 0. We imply below that luq ^ ujq for A < 0. 

One leads to the following results after the corresponding calculations. Two regimes should be considered in this 
case: T <C SCjo and T SCjq. We have at T <C Su^o for RerJ(w,k) expression (pO| which should be multiphed by 
(Do/wq. As a result the spin- wave gap has the form 



D 



(C4) 
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i.e., one leads to Eq. ([3T|) multiplied by y^wo/wo. At T » Sojq we obtain Eq. ([32]) in which the last two terms should 
be multiplied by wq/'^o- In the case of large S and T discussed in Sec. |V]we have for the spin- wave gap expression 
([57)1 multiplied by -^/cjo/i^o- 
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Thermal corrections to the gap proportional to ijj'QT\n{T / {Sujo)) / D'^ from the Hartree-Fock diagram and from the loop 
diagram discussed in Sec. IHI Bl cancel each other. As a result thermal correction to the square of the gap for S ~ 1 is 
much smaller, loIt^''^ / D^"''^ , (see discussion in Sec. IHI C|) . To find it one should take into account O(k^) corrections in the 
expansion of the dipolar tensor Q^^ which depend on the lattice type. 

Although Tc given by Eq. (|5H) can differ several times from the real value of the Curie temperature,^^ this precision is 
sufficient for our estimations. 

Notice that even at large enough dipolar characteristic energy ujo (remaining small compared to the exchange value) , say at 
ujo « O.IJ, we have S > In (^4TrS{D/[SLuo]f^^^ only for pretty large spins S > 8. 
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a) b) c) 




FIG. 2: Diagrams of the first order in 1/S for self-energy parts discussed in this paper. Diagrams (a) and (c) stem from 
three- magnon terms (|10|) in the Hamiltonian whereas (b) comes from four-magnon terms 




k 

FIG. 3: Renormalized spin-wave spectrum at </ik = 0, = 7r/4 and = 7r/2 for a 2D FM at T ^ Sujo with uio = 0.05J and 
S — 1/2. The bare spectrum ai (j)]^ — (solid line) is also presented for comparison. 




FIG. 4; The region in k-plane is shown in which the spin- wave damping is zero at T = 0. This region is defined by Eq. (|47|l . 
The linear dimension of the region is of the order of [Suoo/Df''^. The length of the region along z axis at fc^ = is of the 
order of {Sloo/ D)'^^'^ . As is shown in Fig.[T] z axis is directed along magnetization. Great thermal enhancement of the damping 
inside this region at T S> Suoo is illustrated by Figs. [5] and |6] 
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FIG. 5; Sketch of the spin-wave damping Fk versus the momentum at fc ^ 1 for | sin20kl ~ 1 (solid line), | sin<^k| = 1 (dashed 
line), and | sin0k| = (dotted line). The corresponding dependences of Fk on T, luq and k are also indicated. We imply S ~ 1, 
and T ^ Sujq. Curves for | sin</)k| = 1 and | sin^kl = differ slightly only in the interval A/x/TD < fc ^ Sujq/D. Notice that 
the ratio Fk/tk rises with decreasing k for any given 0k at fc > lS./\/TD, there is the peak at fc ~ A/\/r_D and the exponential 
decay at fc < IS.I\jTD. The ratio Fk/tk is shown in the vicinity of the peak in Fig. (6] 




1 2 3 4 5 

K 

FIG. 6: (Color online.) The ratio of the spin- wave damping and the real part of the spectrum Fk/ek versus the reduced 
wave- vector /t at sin<j!>k = for the classical and quantum 2D FMs on the simple square lattice. The ratios Fk/ek are given 
by Eqs. (|63|l and (|53p in which k = k^ p / [C-y^w'^T) and n = k^J SDcooce/A^ for the classical and quantum FMs, respectively, 
where A> is given by Eq. ((55}. The curve for classical spins is for T <^ j and w <^ j (see Eqs. (|61|l ). There are peaks at 
K ~ ^JSuJoot/T in quantum 2D FMs. Inset: the same is shown for S = 1/2 and S = 3 on a large scale. For all quantum 
magnets the boundary of the area shown in Fig. [3] in which the spin- wave damping is zero at T = is located at k > 10. 



21 



0,22 




0,15 I ■ 1 ■ 1 ■ 1 ■ 1 

0,0 0,5 1,0 1,5 2,0 
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FIG. 7: The peak height on the curve for Fk/fk given by Eq. (|59p versus the dimensionless spin-wave stiffness D/{SJ) for 2D 
FM with S ^ 1 and with the next-nearest-neighbor exchange coupling J' (see the inset). One has D = SJ(1 + 2 J'/ J) in this 
case and it is implied that D ^ Su>o and 8S{J + J') <^ T <^ Tc oc SD even for small D, where 8S{J + J') is the spin-wave 
band width. 
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FIG. 8: The value S ~ (S^) given by Eq. (|66|l versus the dimensionless temperature T/D for the square 2D FM with ujq ~ O.IJ. 




FIG. 9: Some diagrams of the second order in 1/S. 



